Chimera order in spin systems 
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Homogeneous populations of oscillators have recently been shown to exhibit stable coexistence 
of coherent and incoherent regions. Generalizing the concept of chimera states to the context of 
order-disorder transition in systems at thermal equilibrium, we show analytically that such complex 
ordering can appear in a system of Ising spins, possibly the simplest physical system exhibiting 
this phenomenon. We also show numerically the existence of chimera ordering in 3-dimensional 
spin systems that model layered magnetic materials, suggesting possible means of experimentally 
observing such states. 
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Transition to states characterized by simple or complex 
ordered patterns is a phenomenon of central importance 
in equilibrium statistical physics as well as in dynamical 
systems far from equilibrium [l|, Q ■ Examples of simple 
ordering at thermal equilibrium include the aligned ori- 
entation of spins in Ising-like systems, while, in the con- 
text of nonlinear dynamics, this may be observed in the 
phase synchronization of coupled oscillators. However, 
more complex ordering behavior may also occur in vari- 
ous systems under different conditions, especially in the 
presence of heterogeneities. A surprising recent finding is 
that even homogeneous dynamical systems can exhibit a 
robust, partially ordered state characterized by the coex- 
istence of incoherent, desynchronized domains with co- 
herent, phase locked domains Such chimera states 
have so far been observed only in different types of oscil- 
lator populations, including complex Ginzburg-Landau 
equations, phase oscillators, relaxation oscillators, etc., 
arranged in various connection topologies 0-0] ■ Given 
that "chimera" refers to the co-occurrence of incongruous 
elements, one can generalize the concept of chimera-like 
states to include those characterized by simultaneous ex- 
istence of ordered and disordered regions in an otherwise 
homogeneous system. If such a state can occur as the 
global energy minimum of a system in thermal equilib- 
rium, it may widen the scope of experimentally observing 
chimera-like order in physical situations. 

It is with this aim in mind that we investigate chimera- 
like ordering in systems at thermal equilibrium. Specif- 
ically, we consider spin-models as they are paradig- 
matic for different complex systems comprising interact- 
ing components which can be in any of multiple discrete 
states. For example, simple Ising-like models consisting 
of binary-state elements are versatile enough to be used 
for understanding processes operatingin a wide range 
of physical (e.g., magnetic materials |8l4l0|b biological 
(e.g., neural networks [ll|) and social (e.g., opinion for- 
mation [13, EH) systems. The nature and connection 
topology of the interactions between the spins decide 
whether the entire population reaches a consensus, corre- 



sponding to an ordered state, or is in a disordered state 
that corresponds to the stable coexistence of contrary 
orientations. The existence of a chimera state in such 
situations would imply that even though every spin is in 
an identical environment, different regions of the system 
will exhibit widely different degrees of ordering. 

In this paper we report for the first time the occur- 
rence of chimera order in spin systems. This is charac- 
terized for a system of Ising spins by the occurrence of 
regions having different ordering behavior as measured 
by the magnitude of magnetization. The specific system 
we consider in detail is globally coupled and comprises 
two sub-populations (or modules) with the nature of in- 
teractions between spins depending on whether they be- 
long to the same or different groups. Our central result 
is that when subjected to a uniform magnetic field at a 
finite temperature, one of the sub-populations can be- 
come ordered while the other remains disordered. This is 
surprising as both the interactions as well as the external 
field for the two modules are identical. Moreover, the 
chimera state is not a metastable state, but rather the 
global minimum of free energy for the system. The crit- 
ical behavior of the system associated with the onset of 
chimera ordering is established in this paper by an exact 
analytical treatment. We also numerically demonstrate 
the existence of similar complex ordering phenomena in 
three-dimensional spin systems with nearest neighbor in- 
teractions. This opens up the possibility of experimen- 
tally observing chimera states in layered magnetic sys- 
tems, e.g., manganites 

We consider a system of 2N globally coupled Ising 
spins arranged into two sub-populations, each having N 
spins, at a constant temperature T and subjected to a 
uniform external magnetic field H(> 0). A dynamical 
system analogous to our model has recently been ana- 
lyzed by Abrams et al [7j where two clusters of identi- 
cal oscillators, each maintaining a fixed phase difference 
with the others, was shown to possess a chimera state. 
The interaction between two spins belonging to the same 
module is ferromagnetic, having strength J (> 0), while 
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that between spins belonging to different modules is anti- 
ferromagnetic with strength — J 7 (< 0). It is obvious that 
in the absence of an external field, the modules will be 
completely ordered in opposite orientations at zero tem- 
perature. As temperature is increased, the magnitude of 
the magnetizations for the two modules will decrease by 
the same amount, eventually becoming zero at a criti- 
cal temperature, T c . In the presence of an external field 
H that favors spins with +ve orientation, the module 
having negative magnetization will be subjected to com- 
petition between (i) the field H which attempts to align 
the spins along the +ve direction and (ii) the antiferro- 
magnetic interaction J' which is trying to do the oppo- 
site. For a suitable choice of the parameters J, J' and a 
strong field H > H (where H is a threshold field), as 
the temperature is increased from zero, the spins in both 
modules initially remain ordered and are oriented in the 
same direction. Beyond a certain critical temperature 
T c i, we observe that one module becomes more disor- 
dered relative to the other. As the temperature increases 
further beyond a second critical temperature T c2 , the two 
modules again attain the same magnetization, which de- 
creases gradually with T [Fig. Q] (a-d)]. The phase tran- 
sitions at T c i and T C 2 are continuous, characterized by 
critical exponents a and j3 which can be derived exactly. 
For H < Hq, the spins in the two modules are oriented 
at T = in opposite directions, although having the 
same magnitude. At any finite temperature below T C 2, 
the module whose spins were initially oriented opposite 
to the direction of the field is seen to be more disordered 
than the other module. The same critical exponents as 
in the case of H > Hq are observed for the transition at 
T c i-, beyond which the magnetization of the two modules 
are same in magnitude and orientation. 

For the system described above, the energy for a given 
configuration of spins is 

E = — J S SiaSja + J ^ Si a Sjfj — H S Si a , (1) 



where i,j = 1,2, ■•■N labels the spins in a particular 
module and a, /3 = 1,2 identifies the two modules. Since 
each spin is connected to every other spin, mean-field 
treatment is exact for our effectively infinite-dimensional 
system. Thus, the total free energy of the system can be 
expressed as: 

F{m\, 7X1%) = —aN(m\ + m|) + bNm,\m2 — 
HN(m 1 +m 2 ) + Nk B T[S(m 1 ) + S(m 2 )}, (2) 

where mi and m-i are the magnetizations (per spin) of 
the modules 1 and 2, S(m) = |[(1 + m)log(l + to) + 
(1 — to) log(l — to)] — log 2 is the entropy term, and a — 
J(N — l)/2, b = J'N are system parameters (ks being 
the Boltzmann constant). 

To find the condition for equilibrium at a temperature 
T, the free energy can be minimized with respect to toi 
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FIG. 1: (a) Variation of magnetization per spin of the two 
modules (mi, 7712) with temperature. The free energy land- 
scape corresponding to chimera order at fcsT/ J — 5 (b) shows 
that there are two energy minima for mi 7^ m,2 (the curves 
are iso-energy contours and darker shades correspond to lower 
energy), whereas outside the range [T c i,Tca] there is only 
one energy minima (mo) on the mi = mi line as is seen 
for fcsT/J = 8 (c). This is seen explicitly in (d) when the 
free energy per spin F is observed along the curve of steepest 
descent from mo (for T c \ < T < T C 2) or along the curve of 
slowest ascent (for T < T c i or T > T C 2). (e) In the chimera 
ordered state, the system switches between the two energy 
minima corresponding to the two modules exchanging their 
magnetization states between m> and m< shown for MC sim- 
ulations with TV = 100. In all cases, a — 1 and b — H — 10. 



and to 2 to obtain: 



- 2&TO4 +bm?-H+ 1 ^- log l±H!± = 0, (3) 
2 1 — TOl 

knT 1 + TOo , , 

-2oTO 2 + bmi - H + -4- log = 0. 4) 

2 1 — TO2 

Solving Eq. [3] in terms of Eq. [H we obtain the one- 
dimensional map, 
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1 , IT kBT 1 1 + X 

2ax + H log 

2 1 — x 



(5) 



whose solutions of the form g 2 (x) = g(g(x)) — x give 
the extrema ml and of the free-energy F (Eq. [2]). 
Numerical solution for the extrema values shows that for 
suitable parameter values and H > Hq, the system has 
two critical temperatures T c \ and T C 2. For temperatures 
lower than T c \ and above T C 2 the only fixed-point of the 
map g 2 is the unstable fixed point, g(x) = x, of Eq. ([5]). 
Thus, this solution corresponds to mi = to 2 = mo (say), 
where the free energy F(mi,m2) has a minimum. The 
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value for mo is obtained from 

{-2a + b)m -H+ —— log = 0. (6) 

2 1 — too 

However, in the temperature range T c \ < T < T c2 , there 
are two types of solutions of g 2 : (i) a stable fixed point 
toi = to 2 = too [obtained from Eq. ((6])] corresponding 
to a saddle point of the free energy function, and (ii) 
the pair of unstable fixed points mi ^ ra-i which form a 
period-2 orbit of Eq. ([5]) corresponding to a minimum of 
the free energy F. As one of (mi , TO2) is higher and the 
other low, we obtain a chimera state where one module 
is disordered (to<) relative to the other module (to>). 
As shown in Fig.Q](a), the chimera state occurs through 
subcritical pitchfork bifurcations as the temperature is 
increased above T cl or decreased below T c2 . For H < Ho, 
the system exhibits chimera ordering for T > and it 
has a single critical temperature at T c i above which the 
magnetizations of the two modules become same. 

By observing the free-energy F (mi, 1712) landscape in 
the range < toi,to,2 < 1, we obtain a clear physical 
picture of the transition to chimera ordering. The homo- 
geneous state toi = to 2 = too is a local extremum (i.e., 
dF/dmi t 2 = 0) for the range of parameters considered 
here. However, its nature changes from an energy min- 
imum to a saddle point as the temperature is increased 
beyond T c \ and again changes back to a minimum when 
temperature exceeds T C 2. This is seen by looking at the 
matrix of the second derivatives of free energy per site 
with respect to m\,m2'. 

"L - (f a) ■ < 7 > 

where A = —2a + fcgT 1 5 . The eigenvalues of this 
matrix arc A+ = A + b along the mi = to.2 line and A_ = 
A— b in the direction perpendicular to it (parallel to mi = 
— m 2 line). Below T c \ and above T C 2 both eigenvalues are 
positive indicating that mo is a minimum. The transition 
to chimera ordering occurs in the range T c \ < T < T C 2 
when the smaller eigenvalue A_ becomes negative while 
the other eigenvalue remains positive, indicating that Too 
is now a saddle point. This gives us an implicit relation 
for T c as the temperature where A_ = 0, which gives 

k B T c = 2(2o + 6)(l -mo). 

Numerical investigation of the landscape indicates that 
this transition is accompanied by the creation of two min- 
ima away from the toi = TO2 line [Fig. [T] (b-d)]. These 
minima are symmetrically placed about the toi = TO2 
line [as F(m\,m2) = F(ni2,mi)] and correspond to the 
two coexisting chimera states C\ : mi = m>,TO2 = m< 
and C2 : mi = to<, to-2 = m>. The two minima are sepa- 
rated by an energy barrier A = F(mo, too) — F(m > , m<) 
which for a finite system can be crossed by thermal en- 
ergy. This switching behavior between the two chimera 




FIG. 2: (a) Phase diagram in the magnetic field (H), temper- 
ature (kBT/ J) and anti-ferromagnetic coupling (b) parameter 
space obtained by numerical minimisation of free energy for 
a = 1 with the region in which chimera ordering occurs being 
indicated. A cross-section along H — T plane for b — 15 is 
shown in (b). The broken line indicates H = 16, for which 
the critical temperatures are shown by dotted lines. 

states has a characteristic time r ~ exp(A//csT) which 
is indeed observed from Monte Carlo (MC) simulations 
[Fig. [T] (e)]. Note that, each minima corresponds to a 
state with coexisting order and disorder, and hence is 
unlike the minima seen in phase-coexistence state of sys- 
tems such as metamagnets, where each of the minima 
represent a certain type of order [14] . 

Fig. [2] shows the region in (H — T — b) parameter space 
where chimera ordering is observed in our system as 
obtained by numerical minimization of the free energy. 
Temperature induced transitions are always continuous 
whose exponents are analytically derived below. To in- 
vestigate the critical behavior of the system around T c \ 
and T c2 , we shall use the order parameters: 

Pi = mi — m 2 and P2 — 2too — (mi + m 2 ). 

For T c i < T < T C 2 where the chimera ordering is ob- 
served, as mentioned earlier the free energy minima are 
at mi and m 2 while mo corresponds to a saddle point. 
The order parameters p\ and P2 are non-zero in this re- 
gion and and zero elsewhere. When p\,p2 are small, we 
solve for them using Eqs. (J3j> , (|4]) and dU) by expressing 
mi and 7712 in terms of pi,P2, and obtain 

pi oc| T - T c | 1/2 and | p 2 |oc| T - T c | . (8) 

Thus, as T — > T^ or T — > T^ 2 , the order parameters 
vanish continuously with exponents P — 1/2 for p\ and 
(3 = 1 for p 2 - Similar calculations for the field induced 
transition at finite temperature yield identical critical ex- 
ponents. Note that, at zero temperature the field induced 
transition is of first order and its discontinuous nature 
can be shown exactly by analyzing the free energy. The 
values of the exponents for all continuous transitions have 
been confirmed by us numerically. 
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FIG. 3: (a) Schematic diagram of the 3-dimensional layered 
spin system with ferromagnetic (anti-ferromagnetic) interac- 
tions within layers, J (between layers, J') indicated by con- 
tinuous (broken) lines. In the chimera state alternate layers 
show order and disorder, (b) The time-evolution in MC steps 
of the different layers of a system with 32 planes each having 
128 x 128 spins showing chimera ordering for fcgT/ J — 3. (c) 
The magnetizations of different layers of the 128 x 128 x 32 
spin system at different temperatures. Chimera ordering is 
manifested as different values of \m\ for alternate layers (e.g., 
at k B T/J = 3). In all cases J = J' = 1 and H = 1.8. 



We have also analyzed the critical behavior of the 
specific heat C = -T gp where Fq is the equilibrium 
free energy at a given a, b, H and T. Although it in- 
volves both first and second order derivatives of p\ and 
P2, as the most dominant term is d 2 pi/dT 2 , the diver- 
gence at critical temperature is characterized by expo- 
nent a = 3/2 : C oc T - T cliC2 |" 3/2 . 

While the system we have considered so far has the 
advantage of being amenable to exact analytical treat- 
ment, we have also numerically analyzed spin models 
which are closer to real magnetic materials. We have 
performed MC simulation studies of a three-dimensional 
Ising spin model with nearest neighbor interactions hav- 
ing an anisotropic nature. The system emulates a lay- 
ered magnetic system comprising multiple layers of two- 
dimensional spin arrays stacked on top of each other, with 
interactions along a plane being ferromagnetic (J) and 
those between planes anti-ferromagnetic (—J'). Fig. [3] 
shows chimera ordering in such a 3-dimensional spin sys- 
tem with periodic boundary conditions. Similar behav- 
ior was observed in other systems having different sizes, 
parameters and interaction structure, indicating that 
chimera ordering is a robust phenomenon that should be 
possible to observe in an experimental magnetic system. 

In summary, we have shown the existence of a novel 
complex ordering behavior that we term chimera order 
in analogy with the coexistence of coherent and inco- 



herent behavior in dynamical systems. For a system of 
two clusters of Ising spins, where the spins are coupled 
ferromagnetically (anti-ferromagnetically) to all spins in 
the same (other) cluster, subjected to a uniform external 
magnetic field at a given temperature, chimera order- 
ing is manifested as a much higher magnetization in one 
cluster compared to the other. To illustrate the wider 
implication of our result we can use the analogy of two 
communities of individuals who are deciding between a 
pair of competing choices. The interactions of an agent 
with other members of its own community strongly favor 
consensus while that with members of a different commu- 
nity are antagonistic. Thus, given that every individual 
is exposed to the same information or external environ- 
ment, we would expect that unanimity about a particu- 
lar choice in one community will imply the same for the 
contrary choice in the other community. However, the 
occurrence of chimera order suggests that under certain 
conditions, when given the same external stimulus we 
may observe consensus in one community while the other 
is fragmented. The generality of the chimera state as de- 
fined in our paper suggests that it may be experimen- 
tally observed in physical systems. Our demonstration 
of chimera order in a three-dimensional spin system with 
nearest neighbor interactions indicate that a possible ex- 
perimental example can be layered magnetic materials 
(e.g., manganites) having different types of interactions 
between and within layers 0, [l(| ■ Although the present 
paper looks at the case of two competing choices, it is 
possible to extend the analysis to g-state Potts spin dy- 
namics. Given the wider applicability of spin models for 
studying ordering in different contexts, one can consider 
other connection topologies as well as mesoscopic fea- 
tures such as the occurrence of multiple modules (> 2) 
and hierarchical organization. 
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